Convective instability on a crystal surface 
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The distinction between absolute and convective instabilities is well known in the context of hydro- 
dynamics and plasma physics. In this Letter, we examine an epitaxial crystal growth model from 
this point of view and show that a strain-induced step bunching instability can be convective. Using 
analytic stability theory and numerical simulations, we study the response of the crystal surface 
to an inhomogeneous deposition flux that launches impulsive and time-periodic perturbations to a 
uniform array of steps. The results suggest a new approach to morphological patterning. 
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Instabilities dominate the behavior of many nonlinear 
systems. In hydrodynamics 0] and plasma physics 
it is commonplace to distinguish between an absolute and 
a convective instability. In the first case, an initial per- 
turbation spreads more rapidly than it advects and the 
system evolution is insensitive to subsequent perturba- 
tions. In the second case, an initial perturbation advects 
more rapidly than it spreads and the system evolution re- 
mains sensitive to subsequent perturbations. The latter 
characteristic presents an opportunity for control. In this 
Letter, we show that a strain-induced instability that can 
occur during epitaxial growth is convective and suggest 
a novel approach to morphological patterning based on 
the implied control. 

We consider a regular, staircase-like surface composed 
of flat terraces of average width I, The terraces are sepa- 
rated by straight, parallel, atomic height steps with hor- 
izontal positions x n , where the index n grows in the di- 
rection of negative surface slope. We assume that atoms 
impinge on each lattice site at a rate F. This leads to 
a build up of a finite concentration of adatoms on the 
terraces. Adatoms diffuse on terraces and attach to the 
bottom of steps at a rate K . Atoms can also detach from 
steps towards neighboring terraces. These processes lead 
to step motion which can be described by simple equa- 
tions If the flux is large enough, the steps acquire 
a net positive velocity, inducing vertical growth of the 
crystal by one atomic unit after every step has moved a 
distance of one terrace width. 

The equations of motion for the steps are much more 
complicated if there are long-range interactions between 
steps. An example is the growth of a strained film where 
the lattice constant of the deposited material differs from 
the lattice constant of the substrate. The corresponding 
equations of motion p| can be simplified if we assume 
that diffusion is fast. In this limit, the step velocities are 
given by 
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where fi n is the chemical potential associated with adding 
an atom to the solid at the nth step. The /Lt„ are 
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where a reflects the attractive interaction arising from 
the elastic relaxation around each step on a strained 
layer, and (3 reflects the repulsion arising from the in- 
herent stress of each step. 

One solution to this model is uniform step-flow, where 
every step moves with velocity Ft. Under certain circum- 
stances, this steady state becomes unstable and groups 
of steps bunch together |5j. In this paper, we analyze the 
bunching instability from a new point of view. 

The distinction between an absolute and a convective 
instability is most significant for a problem with at least 
one preferred frame of reference. For our problem, the 
lab frame is one such frame. We will also be interested in 
perturbing the step train by supplementing the uniform 
deposition flux with a very narrow beam of atoms that 
can be moved across the surface. This beam is at rest in 
the source frame of reference. 

In the lab frame, the linear stability of uniform step- 
flow motion against a perturbation of the step positions 
Sx n (t) = e exp [i (n£q — tot)] leads to the dispersion rela- 
tion 



Di ab (q,uj) = -i(uj - F£q + F sin £q) 
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Here M is the number of neighbors a step interacts with 
on each side. In a general frame of reference, D(q,u>) — 
Di a b(q,uj + qvf), where Vf is the velocity of the frame 
of reference with respect to the lab frame. Conventional 
stability theory seeks the complex zeroes u>(q) = u)n(q) + 
iu>i(q) of D(q,u>) for given real q. U!i(q) > is a sufficient 
condition for instability of uniform step-flow. Figure [l| 
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shows wj(q) of our model for different values of a. The 
q = mode is marginal for all values of a. When a > 
is small there are two additional marginal modes with 
q = ±q m . All the modes with — q m < q < q m (except for 
q = 0) are unstable. When a is increased, q m increases 
towards tt/£ and the interval of unstable modes becomes 
wider. 

To distinguish between different types of instabil- 
ity, we consider the long-time behavior of the mode 
with wave-number q° that has a zero group velocity: 
du>n(q)/dq\ q=q o = 0. By definition, the system is abso- 
lutely unstable if u>i > and convectively unstable 
if uoi (g°) < 0. In physical terms, this is equivalent to ex- 
amining the time-asymptotic behavior of a disturbance 
launched by an impulse-type (localized in space and time) 
perturbation. When the instability is absolute, the dis- 
turbance spreads in space more rapidly than it advects; 
an observer at any fixed position sees asymptotic growth 
since uii > 0. When the instability is convective, the 
disturbance advects more rapidly than it spreads; an ob- 
server at any fixed position may see transient growth as 
the disturbance advects by, but finds asymptotic decay 
since ivi < 0. 
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FIG. 1. t^i(q) for real q and different values of a. The val- 
ues of the other parameters in units where 1=1 are F = 10, 
K = 6, p = 1 and M = 299. 

Figure [2] is a space-time plot of step trajectories (in 
the lab frame) determined by a numerical solution of the 
equations of motion with an impulsive perturbation ap- 
plied at t = to a single step. In the laboratory, a 
perturbation of this kind can be generated using the nar- 
row beam source mentioned above The perturba- 
tion has no effect on uniform step flow in the portion 
of the surface labeled Region B in Fig. 0. However, in 
Region A, the perturbation creates a disturbance that 
spreads and amplifies in the direction of step flow. v m i n 
and v m ax are the minimal and maximal group velocities 
(in the lab frame) of unstable Fourier modes (for which 
ui(q) > 0). As it happens, v m i n — for this model so 
the disturbance neither spreads out over the entire crys- 
tal nor advects away from the point where the impulse 



was applied. In other words, step bunching as observed 
in the lab frame is "on the border" between absolute and 
convective instability. This fact can be used to test the 
step-bunching model experimentally because it does not 
depend on any of the model parameters. By contrast, the 
transition between absolute and convective instability in 
hydrodynamic systems typically occurs at a single value 
of the control parameter. 
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FIG. 2. Space time plot of a system of 300 steps with peri- 
odic boundary conditions after the application of an impulsive 
perturbation to a single step at x — and t = 0. Each line 
shows the trajectory of a single step in the lab frame. Only a 
small portion of the system is plotted and the step motion is 
amplified. The choice of parameters for this specific system 
in units where £ = 1 is F = 10, K — 6, a = 0.9, P = 1 and 
M = 10. 

We turn next to the response of the growing crystal to 
spatially localized but time-periodic perturbations pro- 
duced by the narrow beam source. Such perturbations 
generate two types of asymptotic behavior which we will 
call switch- on bunching and time-periodic bunching. If 
the source moves with velocity v s (in the lab frame) in the 
interval v m i n < v s < v ma x, the system is absolutely un- 
stable in the source frame and switch-on bunching occurs 
exclusively. For this situation, the step pattern develops 
analogously to the impulsive case (Fig. ||). However, if 
Vs > Vmax or v s < v m i n , the system is convectively unsta- 
ble in the source frame. Switch-on bunching still occurs 
on one portion of the crystal surface, but, in addition, 
time-periodic bunching (which is sensitive to the nature 
of the forcing) can occur on an adjacent portion of the 
surface (Region C in Fig. ||). 

To determine whether time-periodic bunching does oc- 
cur in the regime of convective instability, it is sufficient 
to examine the asymptotic linear response of the step 
system to a spatially localized, time-harmonic source, 

r, , x ( 4:\n£-(v s -Fe)t} 2 \ 

S n (t) = exp l - LL _ iu 8 t , (4) 
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where u) s , a and v s are the source frequency, width and 
velocity in the lab frame. This is called the signaling 
problem |0,0|. In the source frame (n s = n—(v 8 — F£)t/£) 
we find that 



5x ns {t) oc exp(— iu> s t) 



enp(in s £q — a 2 q 2 ) 
c Di ab (q,LO s +v s q) 



dq , (5) 



where _D; a b is the dispersion relation (||) continued to 
the complex q plane and C is a suitable contour in this 
plane. The zeroes of Di a b(q,uj s + v s q) in the q plane 
dominate the integral. Among these, the most important 
zero corresponds to the single mode whose wave-vector 
q*(u s ) has a real part g^(w s ) in the interval [— Tt/£,ir/£] 
and an imaginary part g| (u) s ) that can change sign as uj s 
changes. 

The main result is that there exists a critical frequency 



lo c = \Fsm£q m + q m (v s - F£)\ 



(6) 



If \u> s \ > oj c , the amplitude of time periodic step bunching 
decays as the distance from the source increases. The 
source has little effect on the step-flow pattern in this 
case. However, if \u> s \ < w c , the source induces time- 
periodic step bunching that grows exponentially in space: 



Sx ns (t) oc 
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dDi ab (q,uJ s +v s q) 
dq 



(7) 
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There are two cases to consider. If v s > v ma x, the dis- 
turbance grows in the direction opposite to step flow be- 
cause q*j(uJs) > (Region C of Fig. ||(a)). Conversely, if 
v s < v m i n , we find qj(uj s ) < and the disturbance grows 
along the direction of step flow (Region C in Fig. ^(b)). 
Regions A and B correspond to switch-on bunching and 
uniform step flow similar to the corresponding regions in 
Fig.|. 

For a step-bunch that grows from a time-harmonic per- 
turbation, there is very little nonlinear distortion of the 
bunch shape close to the source, as expected. Frequency 
spectra collected at different spatial locations in Region 
C show that higher harmonics contribute more as the dis- 
tance from the source increases. Nevertheless, the ampli- 
tudes of the fundamental and all higher harmonics satu- 
rate for distances sufficiently far from the source. 

The stabilizing influence of nonlincarity in the step- 
flow case prevents the system from wandering too far 
away from the linear character of the imposed perturba- 
tion. This provides an opportunity to exploit the bunch- 
ing instability to intentionally pattern the crystal surface 
in Region C. The idea is to apply a perturbation prepared 
as a superposition of terms of the form (^) with frequen- 
cies in the range \lu s \ < lu c . As long as nonlinear effects 
can be ignored, each of these terms evolves according to 
Eq. (fjj) in Region C. We can therefore tune the values 
of the amplitudes and phases of the various terms of the 



perturbation, so that at a specific time the step config- 
uration in Region C would be close to a pre-designed 
morphology. 



a) 



250 



200 



150 



100 



b) 



-100 



-150 



-200 




100 



FIG. 3. Space time plots of systems of 300 steps with pe- 
riodic boundaries perturbed by a narrow harmonic source. 
Each line shows the trajectory of a single step in the lab 
frame. Only a small portion of the system is plotted and the 
step motion is amplified. We have indicated the rays which 
correspond to the source velocity and the velocities v m i n and 
v m ax- When the source velocity v aource > v max and \w s \ < uj c , 
periodic step bunching is amplified in Region C in the direc- 
tion opposite to step flow (a). When v aource < v m i n and 
\w B \ < U] c , periodic step bunching is amplified in Region C 
along the direction of step flow (b). The choice of parameters 
for these specific systems in units where I = 1 is F = 10, 
K = 6, a = 0.9 and /3 = 1. The source width is a — I. 

In order to demonstrate the applicability of this idea, 
we attempted to induce a groove-like pattern in a region 
which contained 125 steps. To construct this pattern, we 
used the linear analysis to optimize a small set of ampli- 
tudes and phases for waves with frequencies in the range 
\lo s \ < lu c . We then numerically solved the step equations 
of motion with the designed source. The resulting sur- 
face height as a function of position is shown as circles 
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in Fig. [|. In Fig. |J(a) we have indicated regions analo- 
gous to Regions A, B and C of Fig. |||(a). Figure ||(b) is 
a magnification of the section marked by a two-sided ar- 
row in Fig. |](a). The region of 125 steps we attempted to 
pattern is marked. The solid line in Fig. [|(b) shows the 
predicted linear response of the surface to the designed 
source. Inside the patterned region it is very similar to 
the desired pattern. Near the source (x/l = 1200), the 
shape of the surface obtained from the numerical solution 
of the step equations of motion closely follows the lin- 
ear dynamics. Further from the source, nonlinearity acts 
and we observe a regular sequence of grooves which are 
recognizably "echoes" of the original pattern for many 
periods. We have checked that this behavior is robust in 
the presence of deposition shot noise. 



flow growth. A variety of step bunching scenarios arise 
when conventional step-flow is perturbed by a beam of 
atoms whose flux can be controlled as a function of space 
and time. In particular, there is a regime of time-periodic 
bunching that can be used to launch a sequence of pre- 
designed step-bunch patterns. The nonlinearity of the 
model is such that the bunches do not distort apprecia- 
bly as growth proceeds. 
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FIG. 4. Surface height (in units of the height of a single 
step, ho) as a function of position evolving from a superpo- 
sition of sinusoidal disturbances, h = corresponds to the 
plane x n — nl and the unit o f length is the average step 
seperation on this plane, I = y/£? + hg. The source velocity 
is v s = Ft and it is located at x/l = 1200. Regions A, B 
and C in (a) are analogous to the corresponding regions in 
Fig. ^|(a). (b) is a magnification of the region marked by the 
two-sided arrow in (a). The solid line in (b) shows the pre- 
dicted linear response, while the actual surface morphology, 
resulting from the solutions of the step equations of motion, 
is shown as circles. 
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In a typical experimental system the dispersion rela- 
tion is not known. Nevertheless, one can in principle 
investigate the response of the step system to sources of 
different frequencies experimentally. For each frequency, 
one can measure the induced wavelength (q|j) and am- 
plification rate (gf), as well as the prefactor multiplying 
the exponential in Eq. (^). This information is sufficient 
for the implementation of the design procedure outlined 
above. 

In summary, we have demonstrated the convective na- 
ture of a step-bunching instability that occurs in a re- 
cently proposed model of epitaxial, strain-induced, step- 
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